function ejer1a()

    ptos = load("datos1a.mat");
    puntosX = ptos.xd;
    puntosY = ptos.yd;
	
	plot(puntosX, puntosY, 'r')
	
	y_ajustado = [];

    [m b] = minimos_cuadrados(puntosX,puntosY);
    for i =1:length(puntosX);
        y_ajustado(i) = m* puntosX(i) + b;
    endfor
    
    hold on
    plot(puntosX, y_ajustado,'b')

endfunction
